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The mass of the top quark M top is interesting both as a fundamental parameter of the standard 
model and as an important input to precision electroweak tests. The Collider Detector at 
Fermilab (CDF) has a robust program of top quark mass analyses, including the most precise 
single measurement, M top = 173.4 ± 2.8 GeV/c 2 , using 680 pb _1 of pp collision data. A 
combination of current results from CDF gives M top = 172.0 ± 2.7 GeV/c 2 , surpassing the 
stated goal of 3 GeV/c 2 precision using 2 fb _1 of data. Finally, a combination with current 
DO results gives a world average top quark mass of 172.5 ± 2.3 GeV/c 2 . 



1 Introduction 



The top quark, as the heaviest known fundamental particle, is especially interesting because of 
its large mass. Top quark loops contribute large radiative corrections to other observables such 
as the W boson mass, and the value of those corrections depends strongly on the top quark 
mass. In particular, precise measurements of the top quark and W boson masses are needed 
to constrain the mass of the putative Higgs boson, and for consistency studies if the Higgs is 
observed. The Yukawa coupling of the top quark to a standard model Higgs is roughly one, 
possibly indicating a special role for top in electroweak symmetry breaking. Finally, the mass 
of the top quark has inherent interest as a fundamental parameter of the standard model. 



2 Top Quark Mass at CDF 

The Collider Detector at Fermilab (CDF) is a general-purpose detector at the Fermilab Tevatron, 
a pp collider operating at ^fs = 1.96 TeV. Top quark mass analyses at CDF use events where 
top/antitop pairs are produced strongly. Since top quarks decay before hadronizing to a W 
boson and a b quark, the tt event signature depends on the decays of the two W bosons. Here 



I report on measurements using the "lepton + jets" decay topology, where one W decays to 
an electron or muon and a neutrino, and the other W decays to hadrons; and the "dilepton" 
decay topology, where both W bosons decay to ev e or [iv^. Analyses using the challenging 
"all-hadronic" channel are in progress at the time of this conference; no CDF top quark mass 
analyses select t — > Wb — ► rv T b at this time. In some analyses, the sample purity is significantly 
improved by requiring at least one b tag, namely a displaced secondary vertex inside a jet. 

Extracting the top quark mass from these events is challenging for several reasons. First, 
the events are complicated. Due to the splitting and merging of jets, and the presence of initial- 
state and final-state radiation (ISR & FSR), only in about 50% of selected lepton + jets channel 
events do the four leading jets correspond to the four quarks in the ti decay chain. Even if 
the correct set of jets is selected, it is not easy to determine which jet corresponds to which 
quark; the presence of b tags eases this task. In the case of dilepton events, the two neutrinos 
are observed as a single combined "missing energy" vector in the plane transverse to the beam. 
Second, jets dominate ti events as primary and secondary decay products of the top quarks. But 
jet energies are measured with a relatively poor resolution of roughly 85%/^/Et, contributing 
to statistical uncertainties in the top quark mass; more importantly, there are large systematic 
uncertainties on the jet energy scale that induce the dominant systematic uncertainty on the 
top quark mass. ^ Third, the selected event samples have background contamination at the few 
percent to fifty percent level, depending on the cuts. The background contributions are generally 
well understood, but must be treated properly to minimize associated biases and uncertainties. 

CDF has a robust program of top quark mass measurements. Techniques with different 
sensitivity to event information and systematic effects yield compatible results, giving us confi- 
dence in their accuracy. About ten different analysis techniques have been used at CDF on run 

II data; here we will concentrate on recent results, all using 680-750 pb _1 of data. 

The analysis techniques used by CDF to extract the top quark mass can be separated into 
two categories. Template methods choose some kinematic observable from each event, often an 
event-by-event reconstructed top quark mass. Monte Carlo samples with full detector simulation 
are used to build "templates" of the distribution of this variable for signal samples generated 
with various values of the top quark mass, and for background processes. A likelihood fit 
of the data distribution to probability density functions (p.d.f.'s) derived from the signal and 
background templates yields a measurement of the true top quark mass. Due to this two-step 
process, the chosen observable does not need to be an unbiased estimator of the true mass, and 
in fact the template shapes are often quite complicated (see Fig. [21 for an example). Matrix 
element methods, in contrast, use a direct expression of the probability that a particular event 
is observed given a true top quark mass. The likelihood is built up from parton distribution 
functions, matrix elements for signal and background processes, and "transfer functions" that 
connect quarks at the matrix element level to observed jets, accounting for fragmentation effects 
and detector resolution. The measured top quark mass is then the maximum of the product of 
all the event likelihood curves. Since various approximations must be made in the interest of 
computational tractability, the measurement must in general be calibrated using the observed 
behavior on fully simuluated Monte Carlo samples with known top quark mass. Thus both sorts 
of analyses have as an essential input a good Monte Carlo description of the physics processes 
and a good detector simulation. 



3 Top Quark Mass Analyses With 680 750 pb _1 

In the current data sample of 750 pb -1 , our standard selection cuts ^ select 64 dilepton and 
360 lepton + jets ti candidates. In the 680 pb _1 subset of the data where the CDF silicon 
microstrip detector was operating well, allowing secondary vertices to be reconstructed, we find 
27 dilepton and 252 lepton + jets candidate events with at least one b tag. These are the 
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Figure 1: Results from the two-dimensional template method in the lepton + jets channel. Contours are shown 
at various intervals of log-likelihood in the M top -AjES plane; the maximum likelihood point is given by the red 

cross. 



largest tt datasets yet analyzed, and recent measurements of the top quark mass have statistical 
uncertainties significantly improved over those shown at previous conferences. 

3.1 Two-dimensional Template Analysis in the Lepton + Jets Channel 

The jet energy scale systematic uncertainty is the dominant contribution to systematic uncertain- 
ties on the top quark mass. It is difficult to calibrate jet energies at a hadron collider in general 
because resonant decays to jets are impossible to separate from the continuum background, 
leaving no point of reference for jet energy measurements. But in our reasonably high-purity 
sample of ti events in the lepton + jets channel, there exists a resonant decay of a W boson 
to two jets. CDF has developed a technique to use this information to constrain our jet energy 
scale uncertainty, which results in a lower systematic uncertainty on the top quark mass.^ 

This analysis is an extension of the traditional template method. We would like to measure 
or constrain two quantities, the true top quark mass M top and the difference between our nominal 
jet energy calibration and the true scale. This difference, denoted Ajes> is expressed in units of 
the uncertainty in our nominal calibration, a c , so that we have by definition a prior constraint 
of Ajes = ± 1 <t c . Instead of a single reconstructed variable, we use two observables, the 
reconstructed top quark mass m\ eco and the dijet mass my y The former comes from a kinematic 
fit to the full event, and the latter is the invariant mass of the presumed W daughter jets. 

The quantity m\ eco is primarily sensitive to the value of M top , but is also affected by Ajes- 
The W boson mass is known so precisely that it can be treated as fixed, and the distribution 
of ?7ijj is primarily sensitive to Ajes- But mjj also has a dependence on Mt op , since sometimes 
we erroneously assign one of the b jets from top decay as a W daughter jet. Given all the 
dependencies, we fit simultaneously for the two quantities of interest, using the distributions of 
both observables. In addition, we subdivide the sample into four categories of events, determined 
by the number of b tags and the jet transverse energies. Since the subsamples have different 
signal purity and different resolutions in the two observables, treating them separately increases 
the power of the likelihood fit. 

A global fit to the eight distributions (four subsamples times two observables), using the 
above prior constraint on Ajes, yields likelihood contours in the Mt op -AjEs plane as shown in 
Fig. ^ The anticorrelation arises because a given event sample can be consistent with a high 
top quark mass and a low jet energy scale (i.e. jet energies reconstructed too low) or vice versa. 
The distributions of m\ eco and , along with the predicted distributions at the fitted values of 




Figure 2: Results from the two-dimensional template method in the lepton + jets channel. The plots show the 
distributions of m T t ec ° (left), the reconstructed top quark mass, and (right), the reconstructed dijet mass, for 
each of the four categories of events. In each case, the upper left plot is the subsample of events with two b tags, 
upper right is events with one 6 tag and four high-pr jets, lower left is events with one b tag and a relaxed cut on 
the fourth jet, and the lower right is events with zero b tags, but tighter jet pr cuts. 



M top and Ajes, are shown in Fig. |^1 

The fitted value of Ajes is —0.3 ± 0.6 o~ c , which is consistent with our prior understanding 
of jet energies, but directly translates into a 40% reduction in the corresponding systematic 
uncertainty on the top quark mass. In fact, we have converted the dominant systematic error on 
Mt p into a statistical error, which will improve as we accumulate more integrated luminosity. 
However, the remaining systematic uncertainties still get a large contribution from components 
of the jet energy scale uncertainty that cannot be accounted for with a single jet energy scale 
parameter. Examples are the relative uncertainties in our jet energy calibration between W 
daughter jets and b jets, and between jets of different pt- Using the two-dimensional template 
method, the measured top quark mass is M top = 173.4 ± 2.5 (stat.) ± 1.3 (syst.) GeV/c 2 , the 
most precise single measurement of this important parameter. 



3.2 Two- dimensional Matrix Element Analysis in the Lepton + Jets Channel 

We have recently completed a second analysis that simultaneously measures the top quark mass 
and the jet energy scale, this time based on a matrix element technique. The likelihood to 
observe each event, generically represented as a vector of observables x, depends on the true top 
quark mass, the jet energy scale JES (here defined as a simple multiplicative factor on the jet 
energies), and the signal fraction c s : 

P {x; M topj JES, c s ) = c s P t t{x; M top , JES) + (1 - c s )P w+ict (x; A JES ). (1) 

The meat of the method is in the signal probability: 

P ti (x; Mop, JES) = ~J da tl {y; M top )d qi dq 2 f( qi )f(q 2 )W(x, y, JES). (2) 

Here the differential cross-section a t t, which is determined by the tt production and decay matrix 
element and depends on parton-level kinematic quantities y, is connected with the observables 
x by the transfer function W, which incoporates the dependence on the jet energy scale. The 
variables q\ and q 2 correspond to the incoming partons, and / represents the parton distribution 
functions. The final measurement is based on the product of all the event likelihoods; in this 
analysis the extrinsic jet calibration is not used, so all of the jet energy scale information comes 
from the W mass resonance. 

In this method, the full kinematic and dynamical information from each event is used, in 
contrast to the template method, where each event is reduced to two numbers. However, since a 
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Figure 3: Calibration and result from the two-dimensional matrix element method in the lepton + jets channel. 
The left plot shows the difference between the measured and generated M top in large Monte Carlo samples with 
different generated Mt op (x axis) and different imposed values of Ajes (legend). The right plot shows contours 
at various intervals of log-likelihood in the M top -JES plane; the maximum likelihood point is given by the red 

cross. 



leading-order matrix element is used, events with extra radiation are not well described, and the 
sample is restricted to events with exactly four reconstructed jets. Additionally, a requirement 
of at least one b tag is imposed to improve sample purity, leaving 118 events. 

Since several approximations are made in the likelihoods as written above, we calibrate the 
method against Monte Carlo samples as shown in the left plot of Fig. |H1 Within the statistical 
uncertainty of the Monte Carlo samples used, the method is unbiased. The right plot of Fig. |3] 
shows the likelihood contours, with a result very similar to the template analysis above. Although 
the matrix element analysis uses fewer events, the more powerful method yields a result of 
comparable sensitivity. The systematic uncertainties are dominated by uncertainties in the 
modeling of signal events, for example in the level of ISR and FSR in tt events. Using the two- 
dimensional matrix element method, the measured top quark mass is Mt op = 174.1 ±2.5 (stat.)± 
1.4 (syst.) GeV/c 2 . 



3.3 B Hadron Decay Length Analysis in the Lepton + Jets Channel 

In this template-based analysis we take a different approach to reducing the jet energy scale 
systematic, namely we eliminate it by disregarding the jet energies altogether. Recall that the 
b tagging algorithm reconstructs a secondary vertex inside a jet, identifying the displacement 
between the primary and secondary vertices with the flight path of a B hadron. In ti events, 
the precise magnitude of that displacement depends on the boost imparted to the associated b 
quark, which in turn depends on the mass of its predecessor, the decaying top quark. Thus we 
choose as our observable the decay length in 6-tagged jets; the uncertainties on this quantity are 
uncorrelated with the jet energy scale uncertainty. 

In practice, this is a challenging measurement that is still limited by the statistics of the event 
sample. A proof-of-principle analysis has been performed using lepton + jets events, but dilepton 
and even all-hadronic events could be used in the future. Figure |I] shows the distribution of the 
transverse decay length in the data, compared with the expected contributions from signal and 
background events. The fitted top quark mass is M top = 183.9 (stat.)±5.6 (syst.) GeV/c 2 , 
consistent with other measurements within the large uncertainties. An important point is that 
the statistical and systematic uncertainties are almost completely uncorrelated with other mea- 
surements, so that the method will soon begin to contribute noticeably to M top combinations. 
Studies show that at the LHC, similar methods can make a significant contribution to the overall 
top quark mass uncertainty.^ 
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Figure 4: The transverse decay length L xy , measured in b-tagged jets in the lepton + jets channel, is shown for 
the data compared with the expected contributions from signal and various background processes. 



3.4 Matrix Element Analysis in the Dilepton Channel 

A matrix element technique is also applied to events in the dilepton channel. In dilepton events, 
the top quark mass is harder to reconstruct: the two neutrinos result in an underconstrained 
kinematic system. The dilepton sample is thus particularly amenable to a matrix element 
approach, since the observed quantities and their correlations are used in a more optimal fashion. 

We use a likelihood similar in form to that given above in Section f3. 21 In this case, however, 
the approximations have a more significant effect. The left plot of Fig. El shows the calibration 
curve derived from Monte Carlo samples generated with various values of M% p. The slope of this 
curve is 0.85 ± 0.01, indicating a substantial nonlinearity in the raw measurement, which comes 
primarily from limitations in the background modeling. Additionally, the pull widths for the raw 
measurement are about 1.5 (with no dependence on true top mass), indicating a significantly 
underestimated uncertainty. This effect has been studied extensively, and is understood to result 
from the assumptions that i) lepton energies are perfectly measured; ii) lepton and jet angles 
are perfectly measured; iii) the leading two jets in a signal event always correspond to the two 
b jets from top decay; and iv) the background events are well modeled. After calibration, the 
measurement is unbiased and the error is accurately estimated over a large range of Mt op . 

The right plot of Fig.EJshows the resulting likelihood for the 64 events in the dilepton sample. 
From these events we measure M top = 164.5 ± 4.5 (stat.) ± 3.1 (syst.) GeV/c 2 , the most precise 
single measurement in this challenging channel. The measurement is also performed on the 
subset of 27 events with at least one b tag, yielding M top = 162.7±4.6 (stat.)±3.0 (syst.) GeV/c 2 . 
In both cases, the systematic uncertainty is dominated by the jet energy scale, which contributes 
about 2.5 GeV/c 2 . 



4 Combining Top Quark Mass Measurements 

The most precise information about the top quark mass comes from combining the measurements 
from different analyses, taking advantage of the fact that each method is sensitive to different 
parts of the total available experimental information. We use the BLUE (best linear unbiased 
estimator) method^ to correctly account for the correlations among the systematic uncertainties 
in different analyses. To date, we combine only one analysis from each channel, although work 
is ongoing to understand the statistical correlation between analyses using the same events so 
that more measurements can be included in the combination. 
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Figure 5: Calibration and result from the matrix element method in the dilepton channel. The left plot shows 
the calibration curve of measured vs generated M top in large Monte Carlo samples. The right plot shows the joint 

likelihood curve for the 64 selected data events. 

4-1 CDF Combination 

For this conference, we performed a combination of the CDF run I measurements and the 
measurements described in Sections 13.11 and 111 41 of this document. The combined top quark 
mass is M top = 172.0 ± 1.6 (stat.) ± 2.2 (syst.) GeV/c 2 = 172.0 ± 2.7 GeV/c 2 , exceeding the 
precision of the previous world average. The inputs are consistent when all correlations are 
taken into account: x 2 /ndof = 5.1/4, with a p-value of 28%. CDF's goal for top quark mass 
precision in run Ha (2 fb _1 ) was 3 GeV/c 2 . 

Figure H3 left side, shows an extrapolation of CDF's past performance on top quark mass 
measurements into the future. Two scenarios are given, one where systematic uncertainties are 
fixed and only statistical uncertainties improve with more data; the second where systematic 
uncertainties continue to improve at the same rate as statistical uncertainties, as has been the 
case in the past. The real future of CDF top quark mass measurements probably lies somewhere 
between the two lines. 

4-2 Tevatron Combination 

A combination of current CDF and DO measurements was performed for this conference using 
the same method. © The new world average top quark mass using preliminary measurements 
from both collaborations is M top = 172.5±1.3 (stat.) ±1.9 (syst.) GeV/c 2 = 172.5±2.3 GeV/c 2 . 
The effect on the standard model electroweak precision fit is shown in the right plot of Fig. H3 
which uses the new combination. 

5 Conclusion 

CDF has a strong program of top quark mass measurements, which has recently been augmented 
by a new set of results using 680-750 pb" 1 . A new combination of current CDF results yields 
a top quark mass measurement of 172.0 ± 2.7 GeV/c 2 , putting 1% precision on this important 
parameter in reach as the experiment continues to accumulate integrated luminosity. 




Figure 6: In the left plot, the precision of past CDF top quark mass measurements is plotted as a function of the 
integrated luminosity on a log-log scale. An extrapolation to the future assuming no improvement in systematic 
uncertainties is shown with a solid line; if systematics continue to be reduced as in the past, the dashed line would 
result. The right plot summarizes the current status of the global electroweak fits in the context of the standard 
model, using the new world average top quark mass value produced for this conference. Direct and indirect 
measurements of the top quark and W boson masses are shown using 68% CL contours, and bands indicate the 
preferred Higgs boson mass at each point in the Mtop-Mw plane. 
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